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T’>e  optimlration  of  large-scale  fiyna.nic 
systems  represents  a central  area  of  research 
whose  sijepe,-.:-.ful  outcome  coulh  make  important  con- 
tributions to  the  analysis  of  crucial  national  and 
\orld  problems.  Although  a great  number  of  papers 
have  been  published  on  the  tiieory  of  solving  large- 
scnle  systems,  not  much  in  softvjare  exists  that 
can  successfully  solve  such  systems.  We  believe 
that  there  has  been  little  progress,  because  there 
has  been  little  In  the  way  of  extensive  experimen- 
tation comparing  methods  under  laboratory-like  con- 
ditions on  representative  models.  At  Stanford’s 
Systems  Optimisation  Laboratory  (SOLl,  to  bridge 
this  gup  between  theory  and  application,  we 

(1)  develop  experimental  software  for  solving 
large-scale  dynamic  systems, 

(.?)  systematically  compare  proposed  techniques  on 
representative  models, 

(5)  record  and  disseminate  information  regarding 
experimental  results. 


PILOT  Energy-Economic  Model 

Dynamic  models  that  describe  the  inter- 
actions between  the  energy  sector  and  the  general 
econotsy  help  In  providing  a focus  to  our  research 
in  experimenting  with  large-scale  optimization 
models.  Models  of  this  type  are  under  development 
by  a number  of  groups  to  study  the  energy  crisis 
and  a probable  future  crisis  in  the  area  of  food 
(agriculture).  Developers  of  these  models  could 
make  effective  use  of  techniques  for  solving 
large-scale  systems,  IT  they  were  available. 

Today's  policy  makers  at  industrial, 
governmental  and  international  levels  are  faced 
with  the  decisions  on  providing  the  needed  energy 
ill  th'.-  years  to  come  at  acceptable  social  cost. 
Such  derisions  must  take  Into  account  many  com- 
plex interactions  related  to  the  technology  of 
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energy  .supply,  environmental  side  effects,  energy 
resource  conservation,  etc.,  as  well  as  the  national 
welfare  considerations  of  unemployment,  Inflation 
and  living  standards. 

Some  of  the  important  questions  that  must 
be  considered  in  detail  in  the  formulation  of  the 
energy  policy  (or  policies)  are  the  following: 

(1)  Are  we  using  up  our  cheap  energy  resource.s  too 
quickly? 

(2)  Are  we  making  sufficient  investments  now  so 
that  new  energy  technologies  will  come  into 
commercial  operation  when  needed  in  future 
year.s? 

(3)  Do  we  have  sufficient  physical  capacity  to 
build  the  required  new  plant  and  equipment  in 
the  energy  and  nonenergy  sectors  without 
seriously  hampering  growth  in  consumer  con- 
sumption, or  will  some  sacrifice  In  consump- 
tion be  necessary? 

(U)  VIhat  are  the  various  energy  options  under 
different  patterns  of  crude  oil  Import  price 
realizations? 

(5)  What  will  the  short  and  long  term  Impact  be 
If  oil  and  gas  discoveries  are  less  thaui  pre- 
dicted? 

(6)  Can  we  find  an  energy  policy  that  is  robust, 
l.e.,  one  which  hedges  against  various  con- 
tingencies? 

It  Is  our  belief  that  dynamic  mathematical  pro- 
gramming models  can,  at  the  very  least,  provide 
analysis  and  Information  on  these  and  other  ques- 
tions and  can  substantially  improve  the  understand- 
ing of  the  interactions  that  must  be  considered. 

Such  models  have  been  developed  at  IIASA,  at  the 
Electricity  de  France  and  by  various  groups  in  the 
United  States. 


In  the  Systems  Optimization  Laboratory  at 
Stanford  we  have  under  development  a linear  pro- 
gramming model  for  assessing  energy-economic  op- 
tions in  the  United  States,  called  PIIiOT  (7].  It 
spans  a wide  spectrum  of  activities  of  the  ecotioriy, 
from  exploration  and  extraction  of  raw  energy  to 
industrial  production  and  consumer  demands  for  all 
goods  and  services.  The  data  requirements  there- 
fore cut  across  many  different  sources— consumer 
surveys,  Import/export  and  trade  balance  data, 
manufacturers  surveys,  mining  data,  input/output 
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an  i caplta\  coefficients,  energy  consumption  and 
niil'stitutlon  liata,  enorpy  technology  .lata  fr^m 
Urookhaven  National  Laboratory  and  oil  ami  gas 
expl  ration  and  production  lata.  Hence,  there  Is 
a nontrivial  pstiblom  of  achieving  consistency  and 
of  selecting  a meaningful  level  of  detail  so  that 
the  model  stands  as  a whole  rather  than  as  a con- 
glo'iieratlon  of  parts  that  could  collapse  under 
careful  analysis.  We  believe  that  the  initial  ver- 
sion of  the  model  being  built  will  meet  this  test. 

The  PILOT  model  is  a statement  in  physical 
floir  terms,  to  the  extent  possible,  of  the  broad 
technological  Interactions  within  and  across  the 
sectors  of  the  economy.  Including,  but  in  greater 
detail,  the  energy  sector.  A typical  run  of  the 
model  describes  what  the  country  could  achieve  in 
physical  terms  over  the  long  term,  say  Lo  years, 

A preliminary  version  of  the  model  has 
been  completed  and  several  useful  scenarios  have 
been  run  I IM.  In  1977,  Improved  versions  of  the 
model  will  be  developed,  with  more  detail  regard- 
ing exploration  and  extraction  by  regions,  more 
detailed  modeling  of  various  conversion  processes, 
bettor  rej>resentation  of  foreign  trade,  substitu- 
tion, financial  flows  and  the  effect  of  prices  on 
demand  and  production. 

The  initial  version  of  the  PILOT  model  is 
an  eight  period,  LO  year  model  which  has  approxi- 
mately 00  constraints  and  ?000  variables. 

The  model  Is  a description  in  input/output 
ter"s  of  the  industrial  processes  of  the  economy 
an!  the  demands  for  consumption,  capacity  forma- 
tion, government  services  and  net  exports.  The 
description  of  the  processes  that  provide  useful 
energy  to  the  economy  constitutes  tne  detailed 
energy  submodel.  This  consists  of  technological 
Input/output  deacrlptions  of  the  raw  energy  ex- 
traction and  the  energy  conversion  processes  as 
well  as  the  energy  la?)ort  and  export  activities. 
Four  linkages  Interconnect  the  energy  sector  and 
the  rest  of  the  economy:  energy  demands  of  the 

economy,  bill  of  goods  needed  for  energy  process- 
ing and  capacity  expansion,  total  manpower  avail- 
able to  all  sectors  (Including  energy)  and  a trade 
balance  constraint  which  requires  the  equating  of 
total  exports  to  total  Import.o  when  tliose  Items 
arc  evaluated  In  dollars  over  successive 

five  year  periods.  See  Figure  1,  which  shows  the 
ma.ior  blocks  of  coefficients  in  a time  period, 
and  Its  link  to  the  next  time  period  show  below 
and  to  the  left  of  the  dotted  lines. 

As  noted  earlier,  the  equations  of  the 
model  express  the  balances  of  various  physical 
flows.  For  the  energy  sector,  the  balances  of 
coal,  oil,  gas,  etc.  are  each  expressed  in  BTU 
units.  For  the  economy,  the  units  are  1967 
dollars,  which  are  obtained  by  weighting  the 
underlying  physical  flows  of  goods  and  services 
(assumed  to  be  in  fixed  proportions)  by  1967  prlcaa 
per  unit. 

The  industrial  sectors  of  the  economy  are 
represented  by  a 23-order  input /output  matrix.  The 
sectors  are  grouped  as  follows:  5 energy  sectors, 

1 agriculture,  1 nonenergy  mining,  'j  energy  inten- 
sive manufacturing,  4 energy  nonlntenslve  manufac- 
turing, 1*  services  and  5 capital  formation. 


C',n;unptlon  is  rno.leled  in  terms  of  the  'sn- 
cuniption  vector  of  the  nverot'e  consumer.  This 
sector  doe.-j  not  have  a fixed  bill  of  goods  per 
capita  but  varies  as  a function  of  a parameter 
represent ln<;  tlie  real  consumption  income  attained 
per  capita.  Ba.sed  on  an  analysis  of  historical 
data,  consumption  of  any  Item  as  function  of  aver- 
age con.sumptlon  Income  is  nearly  linear  |2]. 

Capacities  for  each  of  the  IH  nonenergy 
sectors  and  all  of  tiie  energy  processes  are  differ- 
entiated from  one  another.  The  heterogeneous  cap- 
ital equipment  of  the  nonenergy  sectors  Is  depre- 
ciated, whereas  the  energy  facility  capacities  are 
assumed  instead  to  have  undepreciated,  fixed  life- 
times. Construction  lags  are  used  to  specify  the 
time  it  takes  to  build  new  capacity.  These  con- 
struction lags  meiy  be  chosen  individually  for  all 
Ifl  nonenergy  sectors  as  well  as  for  all  energy 
facilities. 

The  exports  are  treated  as  final  demand 
ltcm.s.  The  Import. s nre  considered  in  two  parts, 
competitive  and  noncompetitive.  The  noncompetitive 
Imports  are  for  those  goods  and  services  for  which 
no  domestic  substitutes  exist.  They  are  treated 
as  a part  of  the  technology  of  the  consuming  in- 
dustrial sector.  On  the  other  hand,  competitive 
imports  of  goods  and  services  for  which  domestic 
substitutes  do  exist  are  treated  as  activities  that 
can  augment  the  domestic  production.  Finally,  a 
trade  balance  constraint  ties  together  the  amounts 
of  all  imports  and  exports.  Typically,  we  have 
assumed  over  a five  year  period  that  the  value  of 
total  exports  be  matched  by  that  of  imports. 

The  labor  force  is  assumed  to  be  exogenously 
given.  Also,  average  labor  productivity  is  assumed 
to  grow  at  an  exogenously  given  rate.  In  sample 
runs,  the  "standard  of  attained  appears  to 

be  very  sensitive  to  this  factor, i In  the  base 
case,  2^  per  year  productivity  growth  Is  assumed. 

The  detailed  energy  sector  contains  the  con- 
ventional energy  technologies,  such  as  oil  refinery, 
coal  fired  plant,  etc.,  as  well  as  new  technologies, 
such  as  coal  synthetics,  oil  shale,  plutonlui.i  fired 
reactors,  etc. 

The  description  of  the  energy  sector  includes 
an  accounting  of  reserves  remaining  of  three  ex- 
hsuistible  energy  resources:  oil,  gas  and  uranium. 

For  oil  and  gas,  finding-rate  functions  are  used 
to  specify  the  amount  of  oil  in  place  and  gas  re- 
serves to  be  found  for  a given  amount  of  drilling 
effort.  The  level  of  drilling  effort  Is  endoge- 
nously determined.  The  advanced  (and  expansive) 
techniques  of  secondary  and  tertiary  recovery  are 
also  defined  in  the  model.  Alaskan  oil  production 
and  the  Trans-Alaskan  Pipeline  System  (TAPS)  con- 
struction are  assumed  to  be  exogenously  given.  For 
natural  uranium.  Increasing  facilities  and  man- 
power are  required  to  extract  a ton  of  ore  as  more 
and  more  is  extracted.  In  particular,  progress- 
ively higher  amounts  of  uranium  mining  and  milling 
capacity  are  needed  to  process  the  poorer  grade 
ore  per  pound  of  uranium  oxide  obtained.  While,  in 
principle,  generalized  linear  programming  could  be 
used  to  model  the  nonlinear  functions,  we,  in  fact, 
replaced  the  nonlinear  functions  by  broken  line 
fits. 
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The  I’lLOT  model  can  be  uaed  in  con.junctlon 
with  any  desired  social  ob,1ect,iv.'.  Consideration 
of  reasonable  alternative  ob.lectives  requires 
farther  investigation.  Some  ob.lectives  may  re- 
quire their  expression  partly  in  the  fonn  of  extra 
constraints  as  well  as  in  the  values  of  the  coef- 
ficients of  the  maximand.  The  ob.jective  chosen 
in  the  base  case  in  the  initial  version  of  the 
model  is  the  undiscounted  sum  of  the  gross  na- 
tional consumption  over  nil  I40  years,  subject  to 
(1)  a "monotonic  per  capita"  constraint  which 
states  that  the  average  per  capita  consumption  must 
be  nondecreaslng  over  time,  and  (p)  an  Initial  con- 
dition which  sets  a lower  limit  on  first  period 
consumption.  Experimentation  with  other  maxlmands 
is  possible.  For  example,  we  have  experimented 
with  discounted  gross  national  consumption. 

The  objective  of  PILOT  is  designed  to  per- 
mit one  to  determine  feasible  solutions  to  our 
economy — in  particular,  what  level  of  Investment 
1 in  physical  terms  1 both  in  the  energy  and  non- 
energy sectors  Is  necessary  in  order  to  have  as 
high  a standard  of  living  as  possible  for  the  grow- 
ing population. 

Once  the  physical  flows  are  determined,  it 
is  possible  to  solve  a related  financial  invest- 
ment model.  The  financial  flow  model  calculates 
a system  of  prices,  taxes,  salaries,  profits,  in- 
terest rates,  etc.  that  is  Internally  consistent 
in  the  sense  that  all  economic  agents--consumers, 
producers,  government,  etc. , —receive  sufficient 
monies  to  pay  their  expenses  for  the  specified 
physical  flows.  The  prices  generated  by  the  model 
can  be  adjusted  to  be,  at  the  same  time,  nonin- 
flntlonary,  i.e.,  to  have  the  same  buying  power  as 
base  year  prices.  We  also  are  giving  some  thought 
to  incorporating  in  the  model  production  and  demand 
functions  to  adjust  input/output  coefficients  as  a 
function  of  prices.  In  the  initial  version  of  the 
model  these  coefficients  are  fixed,  however. 

To  illustrate  some  of  the  output  of  the 
model,  a typical  base  case  assumes  a 25t  growth  in 
labor  productivity,  20^  limit  on  the  total  amount 
of  energy  purchased  overseas,  certain  limits  on 
the  rate  of  growth  of  coal  production,  etc.  A base 
case  run  computes  consumption  income  (income  after 
taxes  and  savings).  In  1975  this  income  per  capita 
(in  1975  dollars)  was  about  ^U^oo.  Based  on  these 
assumptions,  the  model  states  it  is  possible  for 
the  country  to  have  a future  consumption  income 
per  capita  relative  to  consumption  per  capita  in 
19^5  as  follows: 

1975  19'0  19^5  1990  19)5  ;000  .005  2010 

1.0  1.0  1.2  1.6.  1.7  1.8  2.0  2.2 

This  possible  future  can  be  compared  with  that  ob- 
tained from  another  scenario,  which  is  the  same  as 
the  base  case  but  restricts  the  use  of  nuclear 
power  plants.  This  naturally  results  in  a lower 
achievable  per  capita  consumption  Income.  Relative 
to  1975,  the  results  are  as  follows: 

kdl  1>80  1'>85  1'>'j0  1995  2000  2005  2010 

1.0  1.0  1.2  1.5  1.5  1.5  1.5  1.5 

Comparing  the  two  scenarios  at  the  yesnr  2010,  we 
have  2.2  vs.  1.5>  i.e.,  the  nuclear  restriction 


could  reduce  the  "standard  of  living"  achieval  i." 
by  OlO  by  JOjt.  This  conclusion  has  been  criti- 
cized becau.se  the  model  assumes  that  consumption 
patterns  of  people  at  different  income  levels  will 
remain  unchanged  (i.e.  they  won't  practice  conser- 
vation or  change  their  life  styles)  and  that  pro- 
duction methods  will  be  no  more  efficient  in  the 
use  of  energy  in  the  future  than  they  are  today. 
This  criticism  we  feel  has  merit  and  we  are,  ac- 
cordingly, considering  revisions  in  the  model  to 
include  more  substitution  and  conservation  possi- 
bilities. 

Because  of  excellent  liaison  with  other 
groups  working  in  the  energy  field,  we  anticipate 
that  the  proposed  physical  flow  mo.iel  will  contri- 
bute to  the  formulation  and  solution  of  the  more 
detailed  specialized  models  under  development 
elsewhere.  In  particular,  the  PILOT  model  is  one 
being  compared  with  other  models  by  the  newly 
formed  U.S.  Energy  Modeling  Forum  in  its  examina- 
tion of  the  feedback  effects  from  the  energy  sec- 
tor on  the  economic  growth. 

Solving  Multi-Time-Period  Models* 

Solving  energy  models  by  commercial  linear 
programming  software  is  proving  to  be  expensive. 
Large-scale  techniques,  such  as  those  under  devel- 
opment at  the  Systems  Optimization  Laboratory,  are 
currently  under  test  to  see  if  they  can  help  solve 
these  models  more  efficiently. 

Conceptually,  the  decomposition  principle 
[ 8]  has  proved  to  be  a natural  approach  to  break- 
ing up  large  systems  and  to  decentralized  decision 
making.  So  far,  computational  experience  has  been 
limited,  but  it  is  known  that  several  devices  can 
be  effectively  combined  with  the  decomposition 
principle  to  accelerate  the  iterative  process. 
Classic  research  along  these  lines  can  be  found  in 
the  work  of  Rosen  [17],  Beale  [3],  Cass  (9]>  Bell 
[U],  Abadie  [1],  Bennett  [5],  as  well  as  in  the 
joint  work  of  Wolfe  and  Dantzig  [8].  Areas  of 
SOL  research  include  (1)  intertemporal  models 
with  staircase  structures,  (2)  the  continuous 
simplex  method  for  linear  control  problems  with 
state-space  constraints,  and  (3)  general  large- 
scale  dynamic  nonlinear  problems. 

Recently,  experiments  have  been  conducted 
at  SOL  comparing  the  Decomposition  Principle  ap- 
proach with  a special  variauit  of  the  simplex  method 
known  as  Generalized  GUB,  These  tests  were  limited 
in  nature  but  indicated  that  Generalized  GUB  is 
superior.  Our  research  on  dynamic  systems  is 
therefore  examining  variants  of  the  simplex  method 
as  well  as  special  methods  for  decoupling  stair- 
case and  block-triangular  systems.  See  Figure  2. 

We  plan  to  compare  these  approaches  with  the  nested 
decomposition  algorithm  of  James  K.  Ho  and  Alan  S. 
Manne  for  staircase  systems  [121. 

Staircase  systems  have  historically  proven 
to  be  very  difficult,  usually  requiring  a dispro- 
portionately large  number  of  simplex  iterations 
to  solve. 


* 

This  section  is  based  on  a sumary  prepared  by 
J.  A.  Tomlin. 
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FIGURE  3.  The  Staircase  Structure  of  the  PILOT  Energy  Model 


J.A.  Tomlin  of  SOL  has  experimented  with  One  of  the  more  promising  methods  that 

a partial  decoupling  of  time  periods  within  a have  been  Investigated  for  reducing  solution  time 

model  by  relaxing  the  Intertemporal  constraints.  for  dynamic  models  Involves  several  modifications 

The  expectation  is  that  such  a relaxation  will  to  the  almplex  method  designed  to  take  advantage 

result  In  a more  easily  solvable  model,  whose  of  the  special  properties  and  behavior  of  such 

solution  can  be  used  as  a starting  point  for  the  andels.  The  essential  property  of  Interest  Is  the 

real  model,  using  some  "gradual"  approach  to  re-  tendency  of  the  same  type  of  activity  to  be  basic 

store  the  Intertemporal  constraints.  As  might  be  over  several  successive  time  periods.  It  there- 

expected,  the  results  of  this  approach  are  quite  fore  seems  desirable  to  Introduce  a profitable 

problem  dependent,  and  sensitive  to  the  degree  and  type  of  activity  In  as  suuiy  time  periods  as  possl- 

the  kind  of  relaxation  employed.  It  appears  that  ble  simultaneously.  To  do  this  M.A.  Saunders  and 

tightly  constrained  econosde  planning  sndels,  of  J.A.  Tosdln  have  explored  variants  of  the  reduced- 

the  type  available  to  us  for  these  experiments,  gradient  method  (a  nonlinear  programming  algorithm 

require  a more  sophisticated  approach.  Other  already  laplamented  by  Murtagh  and  Saunders  In 

types  of  staircase  models  are  being  acquired  to  MINDS  (!<•])  on  these  linear  problems  to  change 

further  test  this  idea.  several  nonbasle  variables  simultaneously.  In 
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contrast  to  the  stamiarJ  simplex  method  which 
changes  only  one  nonbaslc  variable  at  a time.  To 
ensure  that  the  correct  nonbaslc  variables  are 
used,  a "special  pricing"  technique  Is  employed. 
Wlien  the  problem  Is  read  In,  similar  activities 
In  different  time  periods  are  Identified  (from  the 
column  or  variable  names)  and  linked  by  a circular 
list.  Thus  when  an  activity  Is  priced  out  and 
found  to  have  a favorable  gradient,  the  correspond- 
ing vectors  in  other  time  periods  can  be  easily 
found  ani  examined,  and.  If  satisfactory.  Included 
as  candidates  to  be  changed.  It  Is  then  possible 
to  make  a step  which  Introduces  an  activity  In 
several  successive  time  periods  simultaneously. 

Preliminary  experiments  with  the  above 
approach  have  led  to  a reduction  of  20-iQifi  in 
Phase  II  Iterations  when  compared  to  the  standard 
simplex  method  applied  to  the  type  of  economic 
planning  models  referred  to  above.  It  Is  clear 
that  many  tactical  variations  of  the  scheme  need 
to  be  studied. 

If  It  is  advantageous  to  bring  In  an  activ- 
ity simultaneously  In  many  time  periods,  then, 
conversely.  It  should  be  advantageous  to  be  able 
to  also  force  an  unprofitable  activity  to  Its 
lower  bound  In  several  time  periods  simultaneously. 
This  Is  rather  more  difficult,  since  one  cannot 
tell  whether  a whole  group  of  variables  can  reach 
their  bounds  while  maintaining  a feasible  solution 
(at  least,  not  without  incurring  a heavy  computa- 
tional cost).  The  approach  we  have  implemented 
Identifies  groups  of  unprofitable  nonbaslc  activ- 
ities close  to  their  bounds  and  forms  a direction 
vector  scaled  in  such  a way  that  If  one  of  these 
variables  reaches  Its  hound  then  all  of  them  do. 
This  method  has  had  some  success  In  reducing  the 
Iteration  count  when  combined  with  the  "special 
pricing"  described  above.  Again  many  variations 
are  possible,  and  very  considerable  further  ex- 
perimentation Is  required  to  refine  the  methodology 
and  expand  on  the  promising  results  achieved  so 
far. 

VAille  we  have  concentrated  on  means  of 
Improving  the  solution  path  (Iteration  count) 
above,  another  means  of  improving  solution  tech- 
niques for  staircase  models  Is  to  speed  up  each 
step  of  the  simplex  algorithm  by  taking  advantage 
of  the  special  structure  of  the  basis  for  such 
problems.  As  early  as  195**,  OIB.  Dantzlg  [6] 
pointed  out  that  such  problems  exhibit  an  "almost" 
souare  block-triangular  basis  structure  which 
could  be  decomposed  Into  a product  of  a true 
square  block-triangular  matrix  and  another  matrix 
with  only  a few  columns  differing  from  the  unit 
matrix.  The  advantage  of  this  procedure  is  that 
square  block-triangular  matrices  can  themselves 
be  very  efficiently  decomposed  to  give  a very 
sparse  factorization  of  the  basis.  A version  of 
this  method,  employing  modern  factorization  tech- 
niques, has  bean  Impleisented  at  SOL  by  A.P,  Perold 
(a  graduate  student)  and  Is  now  being  tested  on 
problems  of  significant  size.  Karly  Indications 
are  that  this  method  of  handling  the  basis  can  in- 
deed be  more  efficient  than  a direct  treatment 
which  does  not  take  the  staircase  structure  of 
models  Into  account.  Perold* a code  Is  based  on 
sol's  LPMl  linear  prograomilng  code,  as  waa  the 
nested  decomposition  code  by  J.K.  Ho  and  A.S,  Mann* 
[12)  for  the  same  claas  of  probleiss.  This  should 


facilitate  comparison  between  the  specialize': 
simplex  on'i  decomposition  approaches  to  these 
models. 

An  alternative  to  all  of  the  above  numer- 
ical treatments  of  discrete  multl-time-period 
models  is  to  attempt  to  solve  the  underlying  con- 
tinuous time  problem,  which  con  be  thought  of  as  a 
linear  control  problem  with  state-space  constraints. 
G.B.  Dantzlg  and  R.E.  Davis  are  Investigating  a 
"continuous  simplex  method"  for  such  problems.  Al- 
though progress  has  been  made,  much  work  remains 
to  be  done. 
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